Theoretical prediction of Jahn- Teller distortions and orbital ordering in Cs 2 CuCl 2 Br 2 
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With the use of the density function calculations we show that the actual crystal structure of 
Cs2CuCl2Br2 should contain elongated in the ab— plane CuCUBr2 octahedra, in contrast to the 
experimentally observed compression in c— direction. We also predict that the spins on Cu 2+ ions 
should be ferromagnetically ordered in ab— plane, while the exchange interaction along c— direction 
is small and its sign is uncertain. 



PACS numbers: 75.25.-j, 75.30.Kz, 71.27.+a 



I. INTRODUCTION 

The system Cs2Cui(Cl,Br)4 attracts a lot of attention 
mainly because of its unconventional magnetic proper- 
ties. Cs2CuBr4 is the only triangular-lattice antifcrro- 
magnet, which exhibits unusual quantum magnetization 
plateau, [l[ while CS2CUCI4 shows field- induced Bose- 
Einstein condensation of magnons. Thus the investi- 
gation of intermediate compositions of Cs2CuCl4_ a; Br a ; 
may not only shed some light on the magnetic properties 
of both compounds, but may also results in the discovery 
of new phenomena. 

The mixed Cs2CuCl4_ a; Br :l; crystal series was success- 
fully grown at 50 °C with orthorombic crystal structure. 
However, at room temperature in the synthesis process 
the stabilization of a new tetragonal phase for 1 < x < 2 
was observed. [|| 

The crystal structure of tetragonal specimens consists 
of the CuCl2 layers stacking in the c— direction and di- 
vided by the Cs and Br atoms, see Fig.[T] It is important 
to mention that Cu ions in one of the layers placed on a 
top/bottom of the void between CuCU plaquettes of an- 
other layer. According ro Ref. i, the CI and Br ions form 
octahedra surrounding Cu with two short, apical, Cu-Br 
and four long, planar, Cu-Cl bonds. Such kind of the 
compressed octahedra are quite untypical for the Jahn- 
Teller Cu 2+ ions with d 9 electronic configuration for the 
dense crystal structures, and the authors of Ref. H men- 
tioned that the real type of distortions may be hidden by 
the multidomain structures. 

There are a lot of insulating materials, where lig- 
and octahedra surrounding Jahn- Teller active metal ion 
turn out to be elongated: KCrF 4 0,0, Cs 2 AgF 4 @,H|, 
K 2 CuF 4 @,[l3, and others, while there are only few sys- 
tems with the opposite distortion. Generally speaking 
there can be different mechanisms, which stabilize elon- 
gated octahedra. [ll| One of the reasons can be the gain 
in the magnetic energy due to more efficient hoppings 



between half-filled x 2 — y 2 -\ike and oxygen p— orbitals in 
the case of the d 9 configuration. Another one is related 
with the features of the elastic interactions. It can be 
shown that the total energy of elongated octahedra is 
lower than compressed if in the expression for the elas- 
tic energy the terms of higher order than quadratic are 
taken into account. [l2j 

In the present paper, using ab-initio band structure 
calculations, wc found that the crystal structure of 
Cs 2 CuCl 2 Br2 (x = 2), which corresponds to the low- 
est total energy, indeed corresponds to the elongated 
CuCl4Br 2 octahedra. This is in contrast to the observed 
experimental structure, Q but support general tendency 
in the Jahn-Teller distortions to stabilize elongated, not 
compressed octahedra. The more detailed structural 
study should be carried out to confirm (or disprove) the 
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FIG. 1: (color online). The crystal structure of the tetragonal 
phase of Cs2CuCl2Br2. Blue balls are Cu ions, light green, 
brown and purple are CI, Br and Cs ions respectively. Ac- 
cording to Ref. @, the CuCUB^ octahedra are compressed 
along the c-axis. The image was generated using VESTA 
software. B 
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TABLE I: Total energies in meV per formula unit for dif- 
ferent magnetic configurations for the experimental crystal 
structure. [3| The energy of the FM configuration was chosen 
as zero. 





Total energy, meV 


FM 





AFM-A 


3.7 


AFM-C 


2.9 


NM 


8.8 



predicted lattice and orbital ordering. 



II. CALCULATION DETAILS 

The pseudo-potential PWscf code was chosen for the 
calculations. [13| We used ultrasoft pseudo-potentials 
with the nonlinear core correction and the Perdew-Burke- 
Ernzerhof (PBE) version of the exchange-correlation po- 
tential. [lH In order to take into account strong Coulomb 
repulsion on the Cu sites the GGA+U approximation was 
utilized. [l5[ On-site Coulomb repulsion parameter U was 
chosen to be 7.0 eV, while intra-atomic Hund's rule ex- 
change J H = 0.9 eV. [lj] 

The charge density and kinetic energy cut-offs equal 40 
Ry and 200 Ry, respectively. 144 k— points (6x6x4) in a 
full part of the Brillouin zone for the unit cell, consisting 
of 4 formula units (f.u.), were used in the self-consistency 
course. The structural optimization was performed while 
each component of the force were more than 2 mRy/a.u. 
No symmetry operations were used in the course of self- 
consistency. 

The crystal structure was taken for Cs2CuCl2.2Bri.8.[H 
We used the supercell \/2 x y/2 x 2 cell to allow the sim- 
plest types of the orbital ordering. 

III. CALCULATION RESULTS 

The total energy is known to depend on the type of 
magnetic ordering, and this magnetic ordering may have 
an influence on the lattice distortions through the sta- 
bilization of the particular orbital ordering. [171 ] That is 
why it is important to study possible lattice distortions 
together with the analysis of the magnetic interactions. 

The total energies of the following magnetic config- 
uration were calculated: nonmagnetic (NM), ferromag- 
netic (FM), AFM-A type (when all ions in ab— plane are 
ferromagnetically ordered, while the interaction between 
planes is antifcrromagnctic) , and AFM-C type (nearest 
Cu in the ab— plane are AFM ordered, the next nearest 
neighbors in c— direction are ferromagnetically coupled). 

The lowest in energy in the tetragonal crystal struc- 
ture of Ref. H turns out to be the FM configuration (see 
Tab. |T|). The analysis of the occupation matrix shows 
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FIG. 2: (color online). The total and partial DOS for FM 
configuration in the relaxed crystal structure with alternating 
long and short Cu-Cl bond in the ab— plane. Positive (neg- 
ative) values correspond to spin majority (minority) states. 
The Fermi level is set to zero. 



that the single hole in 3c?— shell of Cu 2+ is localized on 
3z 2 — r 2 — orbital in the FM configuration. Such an or- 
bital filling is obviously a result of the local compression 
of CuCl4Br 2 octahedra in the structure of Ref. y. 

The ground state is metallic for any of the investigated 
magnetic configurations. This is in contrast to the fact 
that the undoped materials CS2CUCI4 and Cs2CuBr4 are 
Mott insulators [III [l!| , and the samples of the interme- 
diate compound Cs2CuCl2.4Bri.6 do not have metallic 
shine. 0] 

The compression of the octahedra in CS2CUCI2B1-2 in 
the c— direction prevents the x 2 — z 2 jy 2 — z 2 orbital or- 
der observed in KCuF 3 [IH] and K 2 CuF 4 where 
Cu ions also have one hole in the e g sub-shell. In order 
to allow the same type of the orbital pattern and com- 
pare total energies of different solutions we tetragonally 
distorted CuCUBr^ in the ab— plane and relaxed crystal 
structure with the constrain to keep the same cell vol- 
ume. In addition we also performed the lattice optimiza- 
tion for the initial non-distorted in ab— plane structure, 
since the experimental structure not necessarily corre- 
sponds to the ground state crystal structure in the cho- 
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TABLE II: Distances and total energies per formula unit for 
different magnetic configurations for the relaxed crystal struc- 
tures. The energy of the FM configuration with two long and 
two short Cu-Cl bonds was chosen as zero, it is 557.7 meV 
lower than the FM solution for the tetragonal not relaxed 
crystal structure. 





Cu-Cl dist., A 


Cu-Br dist., A 


Total energy, meV 


FM 


3.02x2/2.25x2 


2.55 x2 





FM 


2.64x4 


2.45 x2 


285.1 


AFM-A 


3.02x2/2.25x2 


2.55 x2 


-0.8 


AFM-C 


2.64x4 


2.44 x2 


280.6 


NM 


2.64x4 


2.44 x2 


286.4 



sen approximation (GGA+U). In effect we obtained two 
crystal structures corresponding to the same magnetic 
order (relaxed initial experimental and relaxed distorted 
in ab-plane structures). The results are summarized in 
Tab. HH where the second, forth, and fifth lines refer to 
structure optimized staring from the lattice presented in 
Ref. (four equal Cu-Cl distances) and the first and the 
third lines correpond to the results for the structure with 
octahedra elongated in ab-planc. 

The lowest total energy corresponds to FM and AFM- 
A type configurations with strongly distorted in the 
ab— plane CuCLj plaquettes, as shown in Fig. [3J The 
energy difference between these two solutions is tiny and 
may depend on the details of the calculations, but both 
lie much lower (~ 558 meV) in energy than the ground 
state magnetic configuration (FM) for the tetragonal not 
optimized crystal structure. Moreover, it is clearly seen 
from Tab. [XT] that only half of this difference can be com- 
pensated by the ionic relaxation which does not change a 
local symmetry of Cu 2+ ions (i.e. that there are 4 equal 
Cu-Cl bonds in ab— plane). 

The second half of the total energy decrease is re- 
lated to the distortions in ab— plane such that they do 
not change the average Cu-Cl bond length, but cre- 
ate checkerboard order of the long and short Cu-Cl 
bonds. The distortions in the ab— plane are accompa- 
nied by a moderate elongation of the CuCLiBr2 octahe- 
dra in c— direction, Cu-Br bond length increases on ~ 0.1 
A, which is compensated by the Cs-Br bond squeezing. 
Thus, instead of compressed in c— direction octahedra, 
the ionic relaxation rather stabilizes the elongated in al- 
ternating directions in ab— plane CuCLiBr2 octahedra. 

Such type of the distortions results in the x 2 — z 2 /y 2 — 
z 2 orbital order, like in KCUF3, see Fig. [3J This orbital 
pattern leads to a strong ferromagnetic super-exchange 
in ab— plane, which stabilizes FM or AFM-A magnetic 
configurations, which agrees with our total energy calcu- 
lations. 

Finally we present the total and partial density of 
states (DOS) plots for the distorted in ab— plane relaxed 
crystal structure, e.g. for FM configuration [Fig. Due 
to a large on-site Coulomb repulsion Cu— 3d band goes 




FIG. 3: (color online). The sketch of the proposed orthorom- 
bic distortions corresponding to the lowest total energy in 
Tab. [HI when all CuCUBr2 octahedra are elongated in alter- 
nating directions in the ab plane (long Cu-Cl bonds are shown 
in red) and the orbital order, which is expected for such kind 
of the distortions (hole orbitals of Cu 2+ are shown). 



away from the Fermi level and is placed mainly in the 
energy range from -5.5 eV to -3 eV. The Br— p and CI— p 
states are concentrated from -3 eV to 1 eV. It is im- 
portant to note, that in contrast to a naive expectation, 
the largest contribution to the bottom of the conduction 
band (as well as to the top of the valence band) comes 
from CI— p and Br— p states: ~3.6 states/ [cV*f.u.], while 
Cu— d provides only 1.9 states/ [eV*f.u.] (the rest be- 
longs to other states of Cu, CI and Br ions). This means 
that the hole is actually localized not on the Cu-x 2 — y 2 
like orbital, but rather on the Wannier orbital, which 
is centered on the Cu site, has x 2 — y 2 symmetry, but 
also has significant contributions (tails) on the surround- 
ing CI and Br ions. This is similar to the situation in 
Cs2Au2Cl6,[20j but here the spatial orientation of the 
Wannier orbital is different. The strong admixture of 
CI— p and Br— p states to the conduction band may be 
the reason of the small band gap value ~0.5 eV in the 
present GGA+U calculations, since these states act as a 
ballast. We expect that the LDA+Uvkf approximation, 
where U is applied not only on the d— part, but on a 
whole Wannier function will result in a larger band gap. 

Summarizing, on the basis of the ab-initio calculations 
we have shown that the layered material Cs2CuCl2Br2, 
containing Jahn- Teller ion Cu 2+ , which was considered 
earlier as a rare example of the Jahn- Teller system with 
localized electrons and e g degeneracy with locally com- 
pressed ligand octahedra, must in fact have elongated oc- 
tahedra with the long axes alternating in the basal plane. 
Thus, yet one more Jahn- Teller material turns out to be 
not an exception, but rather follows the general rule that 
the octahedra around such ions are elongated. This form 
of the Jahn- Teller distortions and orbital ordering should 
lead to strong ferromagnetic exchange in ab-plane, the in- 
terlayer exchange being very weak. The predicted crystal 
and magnetic structure should be observable by the de- 
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tailed structural and magnetic studies. 
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